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We measured N2 fixation rates from oceanic zones that have traditionally been ignored 
as sources of biological N2 fixation; the aphotic, fully oxygenated, nitrate (NOg)-rich, 
waters of the oligotrophic Levantine Basin (LB) and the Gulf of Aqaba (GA). N2 fixation 
rates measured from pelagic aphotic waters to depths up to 720 m, during the mixed 
and stratified periods, ranged from 0.01 nmol N L~ 1 d~ 1 to 0.38 nmol N L~ 1 d~ 1 . 
N2 fixation rates correlated significantly with bacterial productivity and heterotrophic 
diazotrophs were identified from aphotic as well as photic depths. Dissolved free amino 
acid amendments to whole water from the GA enhanced bacterial productivity by 2-3.5 
fold and N2 fixation rates by ~2-fold in samples collected from aphotic depths while in 
amendments to water from photic depths bacterial productivity increased 2-6 fold while 
N2 fixation rates increased by a factor of 2 to 4 illustrating that both BP and heterotrophic 
N2 fixation were carbon limited. Experimental manipulations of aphotic waters from the LB 
demonstrated a significant positive correlation between transparent exopolymeric particle 
(TEP) concentrations and N2 fixation rates. This suggests that sinking organic material 
and high carbon (C): nitrogen (N) micro-environments (such as TEP-based aggregates or 
marine snow) could support high heterotrophic N2 fixation rates in oxygenated surface 
waters and in the aphotic zones. Indeed, our calculations show that aphotic N2 fixation 
accounted for 37 to 75% of the total daily integrated N2 fixation rates at both locations 
in the Mediterranean and Red Seas with rates equal or greater to those measured from 
the photic layers. Moreover, our results indicate that that while N2 fixation may be limited 
in the surface waters, aphotic, pelagic N2 fixation may contribute significantly to new N 
inputs in other oligotrophic basins, yet it is currently not included in regional or global N 
budgets. 
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INTRODUCTION 

In many of the world's oceans, the availability of dissolved inor- 
ganic nitrogen (N) controls primary production in surface waters 
(Falkowski, 1997; Karl et al., 2002). This is because the most 
abundant form of N in the oceans, dinitrogen (N2), is bio- 
logically unavailable to the majority of organisms that reside 
there. Subgroups of prokaryotic organisms (diazotrophs) fix N2 
and convert it to ammonium via the nitrogenase enzyme com- 
plex. This biological process is an important source of "new" 
N to oceanic systems that can stimulate production and growth 
of microbial communities. Current models suggest that new N 
inputs from N2 fixation by cyanobacterial diazotrophs in the 
upper photic layers of the oceans is about 130-140 Tg N y _1 
(Gruber, 2008; Eugster and Gruber, 2012). 

Globally, on geologic timescales, N inputs via N2 fixation 
are thought to be balanced by N losses via denitrification 
(Falkowski, 1997; Deutsch et al., 2007). However, at present, 
rate measurements suggest that N losses from the ocean via 
denitrification (~245-322 Tg N y _1 ) far exceed N inputs via N2 



fixation and terrestrial sources (Codispoti et al., 2001; Gruber and 
Galloway, 2008). The balance between N2 fixation and denitri- 
fication (including annamox) is considered to critically impact 
the oceanic inventory of bioavailable nitrogen and thus primary 
productivity (PP) in the oceans (Gruber and Galloway, 2008). 

Current models generally restrict oceanic N2 fixation to the 
surface, sunlit, photic layer of the oceans where autotrophic 
cyanobacterial diazotrophs occur (Gruber, 2008). However, 
recent work suggests that prokaryotic N2 fixers are more widely 
distributed geographically than previously thought (Moisander 
et al., 2010; Mulholland et al., 2012) and diazotrophs have 
been identified from more diverse pelagic environments includ- 
ing aphotic pelagic waters in oxygen minimum zones (OMZ) 
(Hamersley et al., 2011; Jayakumar et al., 2012). While it is 
known that diazotrophs inhabit shallow benthic environments 
and microbial mats (Capone, 1983; Bebout et al., 1987; Herbert, 
1999), sea grass communities (McGlathery, 2008), coral reefs 
(O'Neil and Capone, 2008), and even hydrothermal vent com- 
munities (Mehta and Baross, 2006), the small areal extent of 
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these habitats is thought to limit their contribution to marine N2 
fixation (Brooks et al., 1971; Carpenter and Capone, 2008). 

For many years, investigations of marine N2 fixation also 
focused on oligotrophic oceanic regions because the presence 
of fixed N was considered inhibiting to active N2 fixation 
(Mulholland and Capone, 1999; Flores and Herrero, 2005; Ward 
et al., 2009). Oligotrophic regions harbor a variety of cyanobac- 
terial diazotrophs including the filamentous cyanobacteria 
Trichodesmium spp., diatom-diazotroph associations (DDA's), 
and unicellular cyanobacteria such as Croccosphaera and group 
A cyanobacteria (UCYN-A) (Bar-Zeev et al., 2008; Carpenter and 
Capone, 2008; Zehr et al., 2008). However, even Trichodesmium, 
the most widely studied marine diazotroph, can fix N2 when 
ambient inorganic N concentrations are on ~5 U.M (Mulholland 
et al., 2001), conditions rarely reached in most oceanic photic 
realms. N2 fixation of marine diazotrophs has been observed in 
nutrient replete cultures (Mulholland et al., 2001) and systems 
with high dissolved inorganic N concentrations (Voss et al., 2004; 
Rees et al., 2009; Mulholland et al., 2012). Data compilation from 
euphotic, mesopelagic, and benthic studies shows that NO^ and 
NH^ concentrations, up to 30 and 200 u,M respectively, can result 
in N2 fixation rates similar to those found in nitrate depleted 
surface waters (reviewed in Knapp, 2012). Moreover, different 
sensitivities to DIN concentrations, exposure time, and compen- 
sation mechanisms all suggest that dissolved N concentrations do 
not control diazotroph distribution and activities in many marine 
environments (Knapp, 2012). 

In addition to photoautotrophic diazotrophs that have been 
previously considered the dominant N2 fixers in marine sys- 
tems (Carpenter and Capone, 2008), heterotrophic diazotrophs 
from a, p, y, 8, and e-proteobacterial phylogenies, as well as 
methanogenic archaea contribute to marine N2 fixation including 
benthic and pelagic realms, mesotrophic and eutrophic environ- 
ments (Chien and Zinder, 1996; Wu et al, 2003; Zehr et al., 2003; 
Halm et al., 2012). Heterotrophic diazotrophs have been identi- 
fied from NO3 -rich deep, aphotic zones of the oceans (Riemann 
et al., 2010), including OMZ (Jayakumar et al., 2012) and from 
near hydrothermal vents (Mehta et al., 2005). In the subtrop- 
ical Atlantic, cyanobacterial nifH genes dominated the photic 
layer (Langlois et al., 2005); while bacterial nifH genes dominated 
the deep aphotic zones (Riemann et al., 2010). In the ultraolig- 
otrophic South Pacific Gyre (SPG) sequence analyses of nifH 
genes and their expressed transcripts revealed a prevalence of het- 
erotrophic diazotrohic y-proteobacteria (Cluster I — Zehr et al., 
2003) and the unicellular photoheterotrophic Group A cyanobac- 
teria (UCYN-A) down to 200 m- the deep chlorophyll maximum 
(Halm etal, 2012). 

While the presence of diazotrophs in oceanic aphotic zones has 
been demonstrated, to our knowledge, only three published stud- 
ies have actually measured rates of N2 fixation in aphotic environ- 
ments; one study was conducted in the Eastern Tropical Pacific 
(Fernandez et al., 201 1), another in the Southern California Bight 
(Hamersley et al, 2011), and the third in a hydrothermal vent 
community (Mehta and Baross, 2006). The first two studies were 
located in OMZ or hypoxic areas. In the Southern California 
Bight aphotic N2 fixation accounted for up to one-third of the 
total areal rates of N2 fixation (55 u,mol N m~ 2 d _1 out of 



150u,mol N nr 2 d _1 ) (Hamersley et al, 2011). While in the 
nutrient-rich OMZ of the Eastern Tropical South Pacific (ETSP), 
N2 fixation rates from sub-surface oxygen deficient waters con- 
tributed about 5 times more N than N2 fixation from the oxic 
photic layer (Fernandez et al., 2011). Moreover, N2 fixation co- 
occurred spatially and temporally with denitrification in the ETSP 
(Fernandez et al., 2011), refuting the paradigm of spatial and/or 
temporal uncoupling between the two processes. 

In oligotrophic environments such as the Levantine Basin (LB) 
(Eastern Mediterranean Sea, EMS) and the Gulf of Aqaba (north- 
ern Red Sea, GA), primary production in surface waters is limited 
at times by nitrogen (N), phosphorus (P), or co-limited by both 
(Paytan et al, 2009; Krom et al, 2010; Tanaka et al., 2011). 
Denitrification in subsurface waters is absent because waters are 
fully oxygenated (Kress and Herut, 2001; Al-Qutob et al, 2002). 
In the EMS, new nitrogen inputs via N2 fixation in the photic zone 
are generally very low, with rates ranging from 0.01 to 0.3 nmol 
NT 1 d -1 . These rates contribute only 1-2% of new produc- 
tion in the LB (Yogev et al, 2011; Berman-Frank and Rahav, 
2012; Rahav et al., 2013a). Although occasional small blooms of 
Trichodesmium and DDA's occur in the GA, N2 fixation rates are 
still relatively low ranging from ~0.05 to 2 nmol NT'd -1 (Post 
et al, 2002; Foster et al, 2009). In the photic layer of the EMS 
and the GA, diazotroph populations include cyanobacteria as well 
as many heterotrophic microbes (Man-Aharonovich et al., 2007; 
Foster et al, 2009; Yogev et al., 2011). 

In this study we measured N2 fixation rates in the aphotic, oxy- 
genated waters of the oligotrophic LB and the GA. Our results 
demonstrate significant N2 fixation rates equal or greater to 
those measured from the photic layers and a prevalence of het- 
erotrophic diazotrophs. This initial data suggests the aphotic oxy- 
genated zones can be considerable sources of new N supporting 
biological production in oligotrophic waters. 

MATERIALS AND METHODS 
SAMPLING LOCATIONS 

Seawater samples were collected from two representative pelagic 
stations that were sampled several times during the mixed and 
stratified periods; one in the GA in the Red Sea and one in the 
LB in the Eastern Mediterranean Sea (Table 1 and below). Water 
samples were collected using 10 L Niskin bottles mounted on a 
rosette equipped with a CTD (Seabird 19 Plus) and fiuorometer 
(Turner designs, Cyclops7) to measure hydrographic properties 
of the water column, and an in situ PAR sensor (QSP-2300, 
Biospherical Instruments Inc) allowing distinction between the 
photic and aphotic zones. 

Gulf of Aqaba (GA) — Water samples from the GA were col- 
lected aboard the R/V Sam Romberg at a station located at the 
northern tip of the Red Sea (Station A 29° 55'N, 34° 28'E) during 
the winter when the water column was well-mixed (March 2010) 
and twice during the stratified summer period (September 2010 
and July 2012). 

Levantine Basin (LB) — Samples from the LB were collected 
during cruises aboard the R/V Mediterranean Explorer and the 
R/V Shikmona. A station located west of the Israeli coastline 
(Station OUT; 32° 57'N, 34° 45'E) was sampled during peri- 
ods when the water column was well mixed (December 2010, 
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Table 1 | Summary of the sampling dates, locations, and depths. 



Region 


Sampling location 


Sampling 


Bottom 


Photic depths 


Aphotic depths 


Hnrt+U n< A 10/ 

ueptn ot u.i /o 








duration 


depth (m) 


sampled (m) 


sampled (m) 


light (m) 


Red Sea 


Gulf of Aqaba 


29° 55N' 


Mar. 2010 


720 


0,20,50,100 


150 


105 




(GA) 


34° 29E' 


Aug. 2010 




0,50, 100 


200,300,500 


120 








July 2012 




10,85 


160, 300, 500, 720 


100 


Med. Sea 


Levantine 


32° 57N' 


Dec. 2010 


1100 


0,20,50,100 


300, 500 


110 




basin (LB) 


34° 45 E' 


Aug. 2011 




5 


250 


140 








Oct. 2011 




5 


300 


130 








Feb. 2012 




5 


300 


110 








Apr.2012 




5 


300 


110 








Mar. 2013 




5, 80 


250 


100 



Additional details are available in the Material and methods section. Med. Sea, Mediterranean Sea. 



February 2012, March 2013) and stratified (August 201 1, October 
2011, April 2012). 

INORGANIC NUTRIENTS 

Duplicate unfiltered water samples were collected in 15-mL 
acid-washed plastic scintillation vials and immediately frozen at 
— 20°C until they were analyzed with a segmented flow Skalar 
SANplus System Instrument as detailed in Kress and Herut 
(2001). The precisions of nitrate+ nitrite (hereafter NO3), ortho- 
phosphate (PO4) and silicic acid (Si(OH)4) measurements were 
0.02, 0.003 and 0.06 u,M, respectively. The limits of detection were 
0.075 |xM, 0.008 |xM and 0.03 |xM for NO3, P0 4 and Si(OH) 4) 
respectively. 

DINITROGEN (N 2 ) FIXATION RATES 

Rates of N2 fixation were measured in field- collected samples 
(4.5 L) using the 15 N2 bubble addition technique described by 
Montoya et al. (1996) and Mulholland et al. (2006) using highly 
enriched 15 N2 gas (99%). The seawater was dispensed into 4.5 L 
Nalgene incubation bottles, 9 mL of 15 N2 was added, and bottles 
were then placed in an incubator with continuously flowing sea- 
water to maintain ambient surface-water temperatures. Triplicate 
bottles were used for each depth sampled. The incubator was 
shaded with neutral density screening to simulate ambient irradi- 
ance conditions at the depths from which samples were collected. 
For samples collected from the aphotic layer, bottles were incu- 
bated in the dark. Both light and dark incubations were 24-48 h in 
duration. While the bubble method may underestimate N2 fixa- 
tion (Mohr et al, 2010; Grofikopf et al, 2012; Wilson et al., 2012), 
particularly during short incubations (hours), our long incu- 
bation times should have minimized this problem (Mulholland 
et al., 2012; Rahav et al., 2013b). See also the discussion regarding 
underestimation of rates due to the bubble vs. enriched seawa- 
ter method (Mohr et al, 2010). To terminate the incubations, 
water was filtered onto pre-combusted 25 mm GF/F filters (450° C 
for 4h). Filters were dried and pelletized in tin disks and ana- 
lyzed on a Europa 20/20 mass spectrometer equipped with an 
automated nitrogen and carbon analyzer preparation module. N2 
fixation was calculated using a mixing model (Orcutt et al., 2001) 
as described by Mulholland et al. (2006) with N solubility factors 
according to Weiss (1970). 



Depth-integrated N2 fixation rates were based on measure- 
ments made at 3 to 6 water column depths spanning both photic 
and aphotic layers. In months where N2 fixation measurements 
were made at only 2 depths, no such calculations were made (i.e., 
Not Available, N.A; Table 2). Depth-integrations included only 
the portion of the water column sampled (range was 150- 720 m 
in depth) and did not extend to the bottom except for one case 
in the GA during July 2012, when N2 fixation rates were mea- 
sured also in bottom waters (Table 1). The percent contribution 
of the aphotic N2 fixation was calculated by difference after sub- 
tracting out the contribution of photic N2 fixation from the total 
integrated water column N2 fixation rate. 

PRIMARY PRODUCTIVITY 

Photosynthetic carbon fixation rates were estimated using Re- 
labeled bicarbonate (Mulholland and Bernhardt, 2005). Water 
samples were placed in triplicate clear 4.5 L polycarbonate 
Nalgene bottles and amended with 99% NaH 13 CC>3 (Sigma) 
to obtain enrichments of about 1% of the ambient dissolved 
inorganic carbon pool. Bottles were incubated under the same 
conditions as for the 15 N2 incubation bottles described above. 
Parallel dark bottles were also incubated for photic depths. As 
for 15 N2, 13 C incubations were terminated by immediately filter- 
ing the entire contents of incubation bottles onto pre-combusted 
25 mm GF/F filters (450°C for 4h). Filters were stored at -20°C 
and then dried and pelletized in tin disks prior to their analy- 
sis on the Europa 20/20 mass spectrometer. Rates of C uptake 
were calculated using a mixing model (as described above) and 
dark carbon fixation was subtracted from the light carbon fixation 
rates for photic incubations. 

BACTERIAL PRODUCTIVITY 

Bacterial production was estimated using the 3 H-leucine 
(Amersham, Specific activity: 160 Ci mmol -1 ) incorporation 
method (Simon et al., 1990). Triplicate 1.7 ml samples were incu- 
bated for 4-8 h at in situ temperatures in the dark. Triplicate, 
trichloroacetic acid (TCA) killed samples (Sigma T9159) served 
as controls. One ml of high 3 H affinity scintillation cocktail 
(ULTIMA-GOLD) was added to samples and they were counted 
using a TRI-CARB 2100 TR, PACKARD scintillation counter. 



www.frontiersin.org 



August 2013 | Volume 4 | Article 227 | 3 



Rahav et al. 



Marine aphotic oxygenic dinitrogen fixation 



Table 2 | Measured physical, chemical, and biological properties of the photic and aphotic layers from the sampling stations in the Levantine 
basin and the Gulf of Aqaba averaged seasonally for when the water column was mixed (winter) and during thermal stratification (spring, 
summer, and fall). 



Location Depth Period Temp. Salinity N0 2 + PO4 BP PP N 2 fixation Contribution of 

(°C) N0 3 [\iM) (m.M) (ngCL" 1 d" 1 | (|xg C L~ 1 d~ 1 ) (nmol N layer to total N 2 

L" 1 d" 1 ) fixation (%) 



Gulf of Aqaba 


Photic 


Mixed 


22- 


-23 


40.5-40.7 


0.09-0.90 


BLD- 


-0.04 


0.4- 


-0.8 


0.1-3.1 


0.05- 


-0.11 


63 






Stratified 


23- 


-27 


40.5-40.8 


0.04-0.06 


BLD- 


-0.04 


1.5- 


-3.8 


0.3-0.6 


0.35- 


-0.63 


25 


Levantine basin 




Mixed 


17- 


-18 


39.1-39.2 


0.06-0.46 


0.02- 


-0.10 


0.4- 


-0.5 


0.7-2.3 


BLD- 


-0.01 


44 






Stratified 


18- 


-28 


39.2-39.5 


BLD-0.20 


BLD- 


-0.03 


0.4- 


-2.5 


0.4-1 .8 


0.01- 


-0.17 


N.A 


Gulf of Aqaba 


Aphotic 


Mixed 


21- 


-22 


40.5-40.7 


0.80-5.61 


0.04- 


-0.36 


0.2- 


-0.4 


0-0.1 


0.02- 


-0.20 


37 






Stratified 


21- 


-23 


40.4-40.7 


0.06-6.24 


0.04- 


-0.36 


0.5- 


-2.4 


0 


0.07- 


-0.38 


75 


Levantine basin 




Mixed 


14- 


-18 


38.9-39.1 


0.20-4.31 


0.10- 


-0.25 


0.8- 


-2.7 


0.3-0.5 


0.01- 


-0.02 


56 






Stratified 


15- 


-18 


38.9-39.1 


0.10-4.52 


0.10- 


-0.20 


4.6- 


-10.7 


N.A 


0.01- 


-0.24 


N.A 



W.A not available, BLD, below detection limit. 



AMIN0-ACIDS ENRICHMENT EXPERIMENT 

We assumed that enriching the water with amino acids would 
promote conditions that were advantageous to heterotrophic 
bacteria in general and also to heterotrophic diazotrophs. We 
assumed that if heterotrophic diazotrophic bacteria were lim- 
ited by energy (carbon), the addition of high-carbon containing 
amino acids would alleviate this limitation allowing these popu- 
lations to fix N2. The idea was to distinguish general bacterial het- 
erotrophic activity (bacterial production rates) from diazotrophic 
ones (N2 fixation rates). The response of the microbial commu- 
nity to additions of an amino acid mixture (Sigma A9906) was 
tested at 4 aphotic (160, 300, 500, 720 m) and two photic (10, 85) 
depths (Tables 1,2) during July 2012 in the GA. Amino acids were 
added to 4.5 L polycarbonate Nalgene bottles filled with whole 
water from a specified depth at a final concentration of 500 nM 
C immediately prior to making 15 N2 additions. 

PYROSEQUENCING 

DNA was extracted from a representative aphotic depth in the 
LB stations. Fifteen L of seawater were filtered onto 0.2 u,m 
Supor® filter (PALL) and stored in — 80° C until analysis in the 
lab. DNA was extracted using a genomic DNA kit (Roche) and 
the extract was purified with a QIAGEN purification kit accord- 
ing to the manufacturer's instructions. Amplification of nifH 
genes from genomic DNA was performed as described elsewhere 
(Zehr and McReynolds, 1989). The nifH sequences retrieved 
from the pyrosequencing analyses [performed by Research and 
Testing Laboratories (MRDNA, Texas, USA)] were trimmed and 
checked for quality using QIIME (Caporaso et al., 2010). The 
cleaned sequences were compared to a nifH database (Zehr et al., 
2003). In the aphotic layer of the LB only 20 sequences with an 
open reading frame were retrieved. These sequences were fur- 
ther confirmed using BLASTX (Altschul et al, 1990) with a cutoff 
of 95%. The ni/H-like sequences were integrated into a nifH 
phylogenetic tree created with MEGA 4.0 software described in 
details by Tamura et al. (2007). nifH gene sequences were sub- 
mitted to GeneBank and were assigned the accession numbers 
KC222033-KC222052. 



DETERMINATION OF TRANSPARENT EX0P0LYMER PARTICLES (TEP) 
CONCENTRATION AND SIZE 

Triplicate water samples (100 ml) were filtered onto a 0.4 u,m 
polycarbonate filters and stained for 2 s with a 0.02% Alcian Blue 
solution. TEP concentrations (ixg GX equivalents L _1 ) were mea- 
sured according to Passow and Alldredge (1995). A conversion 
factor of 0.51 was used to convert from micrograms of Gum 
Xanthan (GX) equivalents to micrograms of carbon (Engel and 
Passow, 2001). 

For TEP size (area - |xm 2 ) measurements, samples (100 ml) 
were filtered gently (~100 mbar) onto 0.4 [im polycarbonate fil- 
ters, stained for 10 s with freshly prepared Alcian Blue solution 
(0.02%), rinsed with 5 ml of double distilled water (DDW), 
placed onto microscope slides and kept at 4°C. The filters were 
analyzed under X40 magnification after 1-2 days and at least 90 
pictures per slide were taken. The pictures were analyzed with 
image J software (http://rsbweb.nih.gov/ij/) to determine the size 
of TEP. 

TEP ENRICHMENT EXPERIMENT 

The association of aphotic diazotrophs with TEP was explored by 
adding 4 different quantities of pure GX (Sigma G1253) to LB 
seawater taken from 250 m (March 2013). Three replicate bottles 
were tested per treatment and per unamended controls contain- 
ing seawater only. The bottles were incubated in the dark at ambi- 
ent temperatures (15°C) for 24 h, after which they were sampled 
to measure TEP concentration and size and rates of N2 fixation 
and BP. 

STATISTICAL ANALYSES 

Data were displayed as means with corresponding standard 
deviations for correlation and regression analyses between vari- 
ables. Comparisons between unamended controls (no addition) 
and the amended microcosms (+Dissolved free amino-acids, 
DFAA) at the end of the experiment, a One-Way ANOVA fol- 
lowed by Tukey's multiple comparison post hoc using R soft- 
ware (www.r-project.org). P < 0.05 was considered statistically 
significant. 
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RESULTS 

PHYSICAL AND CHEMICAL CHARACTERISTICS OF THE APHOTIC LAYER 

The aphotic layer was determined at each station as the depth 
below which irradiance was <0.1% of surface irradiance. The 
depths of the photic layer varied between stations and sea- 
sons, ranging from 105 m to 140 m at both sampling locations 
(Table 1). Temperature and salinity in the aphotic layer of the 
GA station were 21.5 ± 0.5°C and 40.6 ± 0.1, respectively dur- 
ing both the mixed (winter) and stratified (summer) periods. In 
the photic layer of the GA, temperatures were 22 ± 0.5°C during 
winter and 25 ± 2°C during the stratified months (Figures 1A,B; 
Table 2). At the LB station, temperatures and salinities of the 
aphotic layer were significantly lower than those measured for the 
GA station, averaging 16.1 ± 1.1°C and 39.1 ± 0.1, respectively 
(Figures 1C,D; Table 2), while within the photic layer tempera- 
ture and salinities averaged 23 ± 5°C and 39.3 ± 0.2, respectively, 
during the stratified period (Figures 1C,D; Table 2). 

Average dissolved NO3 concentrations in the nitricline and 
sub-nitricline depths of the GA were 3.0 ± 2.1 (jlM with the 
nutricline ranging between 200-400 m during the winter mixing 
(n = 9, Figure 1A; Table 2). NO3 concentrations remained sim- 
ilar (3.0 ± 2.5 u,M) during the stratified summer for the same 
depths although the nutricline deepened to 150-500 m (n = 13, 
Figure IB; Table 2). PO4 (DIP) concentrations varied similarly 
in both periods, averaging 0.2 ± 0.2 u,M (Figures 1A,B; Table 2). 



Surface NO3 and DIP concentrations were low and close to their 
detection limit during both periods (Figures 1A,B; Table 2). 

Unlike the GA, no seasonal fluctuations in NO3 and DIP con- 
centrations were recorded in the aphotic depths of the pelagic 
LB station (Figures 1C,D). At these depths lower NO3 and DIP 
concentrations were 2.2 ±2.1 u,M and 0.1 ± 0.5 u,M, respectively 
[n = 7-9, Figures 1C,D; Table 2). 

N 2 FIXATION RATES 

Photic N2 fixation rates compiled from both stations, at all 
depths and seasons yielded rates usually below 0.1 nmol N L _1 
d _1 excluding rates measured in the GA during summer (0.35- 
0.63 nmol N L _1 d _1 ) (Table 2). Aphotic N2 fixation rates ranged 
from 0.01 to 0.38 nmol N L _1 d _1 at both locations sampled 
(Table 2). Aphotic rates obtained in the GA during the mixed 
period ranged from 0.02 to 0.2 nmol N L _1 d _1 and averaged 
0.1 ±0.1 nmol N L d _1 while during the stratified period 
rates ranged from 0.07 to 0.38 nmol NT 1 d _1 and averaged 
0.30 ± 0. 10 nmol NT'd" 1 (Table 2). At the LB station, aphotic 
N2 fixation was uniformly low during the mixed period; 0.01 ± 
0.01 nmol N L _1 d _1 and increased ~20-fold during the stratified 
period to 0.24 ± 0.09 nmol N L" 1 d" 1 (Table 2). 

To understand the significance of aphotic N2 fixation relative 
to the total areal rates of N2 fixation, we calculated that during 
our sampling times, N2 fixation from aphotic depths sampled 
accounted for 37 to 75% of total areal N2 fixation at the GA sta- 
tion and ~50% of the total areal N2 fixation at the LB station 
during the mixed period (Table 2). 

PHYT0PLANKT0N AND BACTERIAL PRODUCTIVITY RATES 

PP rates in the GA averaged 3.0 ± 0.1 \ig C L _I d _1 and 0.5 ± 
0. 1 (xg C L _1 d _1 for the mixed and stratified periods, respectively 
(n = 3, PAR at surface of 1000-1500 |i,mol quanta nr 2 s _1 ). 
PP at the LB station was generally lower; averaging 2.1 ± 0.2 |i,g 
Cr 1 d _1 and 1.8 u,g C L d _1 for the mixed and stratified 
periods, respectively (n = 3, PAR at surface of 1500-1800 (xmol 
quanta raT 2 s _1 ) (Table 2). Bicarbonate uptake rates were usually 
below detection limits within the aphotic zone (below 105-140 m, 
depending on sampling location) excluding a wintertime mea- 
surement in the LB station (200 m) when C uptake averaged 
0.4±0.1(igCL" 1 d" 1 (Table 2). 

Bacterial productivity (BP) rates in the photic waters from 
both locations ranged from 0.4 to 3.8 u,g C L _1 d _1 , while in the 
aphotic zone BP rates were 0.2 to 10.7 u,g C L _1 d _1 (Table 2). 
Rates of N2 fixation from both the photic and aphotic layers 
were compared to BP and PP (Figure 2). BP correlated signif- 
icantly and positively with N2 fixation rates within the photic 
zone of both basins CR 2 = 0.70, P = 0.001, n = 12), while there 
was no apparent correlation between N2 fixation rates and PP 
(Figure 2A). Moreover, BP and N2 fixation were also positively 
correlated within the aphotic depths (.R 2 = 0.58, P = 0.003, n = 
10) (Figure 2B). 

To further understand what limits N2 fixation by heterotrophic 
diazotrophs in oligotrophic and aphotic systems we added dis- 
solved free amino acids (DFAA) to microcosm experiments con- 
ducted in the oligotrophic stratified period during July 2012 in 
the GA. In water samples collected from the aphotic layer, DFAA 
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FIGURE 1 I Spatial (depth) and temporal changes of physical and 
chemical parameters at the sampling stations in the Gulf of Aqaba 
(A,B) and Levantine basin (CD) during the mixed (A,C) and stratified 
(B,D) periods. 
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additions enhanced BP rates by a factor of 2-3.5 and N2 fixation 
rates increased by a factor of 1.5-2 (ANOVA, P < 0.05) (Table 3). 
Within the photic layer, the response of bacteria to DFAA addi- 
tions was greater than that observed in the aphotic zone; BP 
increased 2- to 6-fold and N2 fixation rates increased by a factor 
of 2 to 4 (ANOVA, P < 0.05) (Table 3). 



TEP AS SOURCE OF CARBON OR MICROHABITAT 

TEP concentrations ranged from 100 to 5000 (xg GX L _1 and 
varied in size; the two-dimensional TEP size ranged from 24 to 
40 |xm 2 with the largest particles measured at TEP concentra- 
tions of ~3500 |xg GX Lr 1 (Figure 3). BP increased almost 7-fold, 
from ~0.6 (xg C L _1 d _1 in the unamended control bottles to 
~4 |xg C L _1 d _1 in incubations amended with the highest TEP 
concentrations (5000 u,g GX L _1 ). N2 fixation was also stimu- 
lated (~57-fold) by TEP additions increasing from 0.06 nmol N 
L _1 d _1 in the controls to 3.46 nmol NT 1 d _1 in the bottles 
amended with ~3500|ig GX L _1 . N2 fixation rates were pos- 
itively correlated with both TEP concentration (_R 2 = 0.49, P = 
0.07, n = 15) and TEP size (R 2 = 0.90, P = 3 x 10~ 5 , n = 15) 
(Figure 3). 

MOLECULAR DIVERSITY OF nifH 

NifH-Vke genes were retrieved from the aphotic layer of the 
LB and were compared with nifH genes from known bacte- 
rial isolates. Twenty sequences were retrieved and all clustered 
together with heterotrophic bacteria and were most similar to 
Pseudomonas stutzeri (Figure 4). All of our sequences were sep- 
arated from the autotrophic clusters. 

DISCUSSION 

Most marine N2 fixation studies neglect the potential contribu- 
tion of biological N2 fixation in deep, aphotic, NGvrich waters 
that comprise the majority of the world's oceans (Karl et al., 2002; 
Zehr and Ward, 2002). This is likely because phototrophic and 
cyanobacterial diazotrophs have an ecological advantage in sur- 
face waters where productivity is often limited by the availability 
of fixed nitrogen yet there is sufficient energy (irradiance) to drive 
both photosynthetic C uptake and the costly process of N2 fix- 
ation (Howarth et al., 1988; Postgate, 1998; Karl et al, 2002). 
Currently these assumptions are being modified as significant N2 
fixation rates are reported from marine areas with high DIN con- 
centrations (review in Knapp, 2012). Here we provide evidence 
that aphotic, nitrogen-rich waters also harbor active diazotrophs 
that contribute substantial new N to the oceanic N inventory, and 
that N2 fixation by these bacterial diazotrophs may be limited by 
the supply of organic C. 

In this study we have documented for the first time significant 
N2 fixation rates from the aphotic, oxidized waters (Figure l)of 
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FIGURE 2 I The relationship between N2 fixation and primary and 
bacterial productivity within the photic (A) and aphotic (B) layers of 
the LB and GA stations. Values are averages and error bars are standard 
deviations of three independent replicates performed for each control and 
treatment incubation. 



Table 3 | The effect of dissolved free amino acid (DFAA) enrichments (+ DFAA) on bacterial productivity (BP) and IM2 fixation during the 
stratified summer period (GA, July 2012) within different depths of the photic and aphotic layers. 

Sampling depth (m) BP ((ig C L~ 1 d~ 1 ) N 2 fixation (nmol N L _1 d _1 ) 

Control + DFAA Control + DFAA 



Photic 
Aphotic 



10 

85 

160 

300 

500 

720 



0.30 ±0.40 
1.15 ± 0.33 
2.34 ±0.69 
2.20 ±0.70 
0.46 ± 0.07 
2.38 ± 0.43 



2.40 ±0.38 
4.17 ±0.68 
4.68±0.14 
5.32 ±0.20 
1.57±0.19 
4.86 ± 0.47 



0.35 ±0.09 
0.30 ±0.11 
0.38 ±0.05 
0.23 ±0.05 
0.07 ±0.03 
0.16±0.10 



0.76 ±0.10 
1.12 ±0.27 
0.70 ±0.20 
0.34 ±0.02 
0.13 ±0.06 
0.31 ±0.06 



Measurements are based on 3-4 replicates for each sampled depth. 
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FIGURE 3 | The relationship between l\l 2 fixation rates of aphotic 
diazotrophs and transparent exopolymer particles (TEP) concentrations 
(A) and size (B) in experimental manipulations of LB waters from 250 m 



depth to which we added 4 concentrations of pure gum xanthan. Values 
are averages and error bars are standard deviations of three independent 
replicates performed for each control and treatment incubation. 
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FIGURE 4 | nifH gene tree of clones retrieved from the aphotic layer in the Levantine Basin. Numbers in parentheses represent the number of clones in 
the library that were identical to the clones in the phylogenetic tree. Bootstrap values above 50% are indicated above the nodes. 



the oligotrophic LB and the GA. The measured N2 fixation rates 
in the aphotic waters of both LB and GA ranged from 0.01 nmol 
NT 1 d _1 to 0.38 nmol N L _1 d _1 . These rates were similar 
to the N2 fixation rates reported from aphotic depths in the 
hypoxic Southern California Bight (Hamersley et al., 2011), yet 
were ~10-100 fold lower than those measured in the Peruvian 
upwelling OMZ (Fernandez et al., 201 1). 



The significance of our measured rates was realized upon inte- 
gration of water column N2 fixation from surface to the deepest 
depth sampled (Table 1). Thus, the contribution of aphotic N2 
fixation accounted for ~50% from the total N2 fixation in the 
LB station, where winter mixing is shallow and the aphotic water 
hardly mixes with the photic layer, and 37 to 75% within the GA 
station (Table 2). 
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The substantial role of aphotic N2 fixation in the LB and 
the GA may be especially important to basin-wide productivity 
given the low rates of N2 fixation in the photic layers of both 
these nutrient-poor environments (Foster et al., 2009; Berman- 
Frank and Rahav, 2012). Within the photic layer, N2 fixation 
was undetectable in a third of the measurements made in the 
Eastern Mediterranean Sea during several studies (m > 250) span- 
ning seasons and years (Ibello et al, 2010; Yogev et al., 2011). In 
contrast, measureable rates ofN2 fixation were observed in >90% 
of samples collected from the aphotic zone (n = 48), regardless of 
season, thermal stability or depth (Table 2), further emphasizing 
the importance of N2 fixation in these realms and the requirement 
for subsequent comprehensive investigations that could constrain 
modeling efforts. 

Our measured aphotic N2 fixation rates may actually repre- 
sent an underestimate of this process as: (1) our sampling did 
not cover the full depth of the water column — thus integration 
may yield higher contribution of aphotic to the total N2 fixation 
(see methods); (2) we used the common gas bubble enrichment 
method (Montoya et al., 1996 and see Materials and Methods) 
which may underestimate N2 fixation rates by a factor of 2-6 
compared to the enriched seawater method (Grofikopf et al., 
2012; Wilson et al., 2012). Early field comparisons from the LB 
photic waters between the two methods showed a 2-3 fold under- 
estimation by the gas bubble technique (Rahav et al., 2013b). 
Assuming our aphotic rates were underestimated by a similar 
ratio of ~50%, our measured aphotic N2 fixation would increase 
accordingly so that maximal aphotic rates would be ~0.8 nmol N 
IT 1 d~ 1 , and with DFAA addition up to 1 .4 nmol N L~ 1 d~ 1 . Yet, 
as the estimated increase by 50% would also be accounted for in 
photic rates of N2 fixation, the % contribution of aphotic to total 
N2 fixation would remain the same. 

Energetic requirements for diazotrophy are high (Postgate, 
1998), thus we initially assumed that phototrophic (mainly 
cyanobacterial) diazotrophs would dominate the photic zone, 
while at aphotic depths, N2 fixation rates would be extremely 
low and accomplished by heterotrophic bacterial populations that 
would be prone to carbon limitation due to the low overlying 
productivity. However, we found that in the LB and GA sta- 
tions heterotrophic diazotrophs were the dominant N2 fixers in 
both photic and aphotic zones. This was confirmed by examining 
the relationship between N2 fixation and BP and PP (Figure 2), 
molecular analysis of nifH genes (Figure 4), and the stimulation 
of N2 fixation after amino acid amendments (Table 3). In the 
stratified photic zone, a significant correlation was found between 
BP andN 2 fixation (R 2 = 0.70, P = 0.001, n = 10) and both were 
stimulated by the addition of DFAA (Table 3). In contrast, no 
correlation was found between N2 fixation and PP (Figure 2A). 
Although ample solar irradiance was available, heterotrophic bac- 
terial diazotrophs appear to dominate in this system relying on 
dissolved and particulate organic carbon (DOC and POC, respec- 
tively) as energy sources (Goldthwait et al., 2005; Azam and 
Malfatti, 2007; Eichinger et al., 2009). 

Within the aphotic layer, N2 fixation was positively correlated 
with BP (R 2 = 0.58, P = 0.003, n = 10), (Figure 2). Moreover, 
a 2-3 fold enhancement of N2 fixation rates and enhanced BP 
were observed following DFAA amendments (Table 3). DFAA 



are labile dissolved organic carbon (DOC) and nitrogen (DON) 
sources for microbes (Middelburg and Nieuwenhuize, 2001) that 
can be directly utilized by bacterial heterotrophs (Davis et al., 
2009). Our data thus suggests that the addition of C and N 
sources stimulated both diazotrophic and non-diazotrophic het- 
erotrophic populations. Moreover, the increase in N source con- 
comitantly with enhanced N2 fixation rates demonstrate that 
the higher concentration of dissolved N sources (either DON or 
bacterially catabolized into DIN) were not sufficient to inhibit 
diazotrophy. These results highlight the significant role of het- 
erotrophic bacteria as important components of the photic and 
aphotic diazotrophic community and suggest that the availabil- 
ity of organic matter may limit bacterial diazotrophy in these 
oligotrophic regimes. 

Previous molecular characterization of diazotrophs that were 
actively expressing nifH genes from the photic zones of the LB 
and GA demonstrated a diverse array of both heterotrophic bac- 
teria and cyanobacterial phylotypes from Clusters I, II, and III 
(Foster et al., 2009; Yogev et al., 2011). In this study, phyloge- 
netic analyses of nifH from the aphotic depths indicated that 
diazotrophs present were within Clusters I and III, phylotypes 
predominantly comprised of heterotrophic bacteria, including a 
and y proteobacteria (Chien and Zinder, 1996) (Figure 4). The 
nifH phylogenetic diversity obtained from the aphotic zone dur- 
ing this study was similar to that recovered from the Southern 
California Bight (Hamersley et al., 2011) and anaerobic bac- 
teria from the Sargasso Sea (Hewson and Fuhrman, 2007) 
(Figure 4). 

Our results suggest that in the aphotic zones of both LB and 
GA heterotrophic bacterial N2 fixation is a considerable source of 
new N. Yet these environments have not heretofore been consid- 
ered as major sources of new oceanic N inputs, likely because they 
are dark, oxygenated, and N03-rich (Figure 1). Heterotrophic 
diazotrophs in aphotic areas may be carbon (energy) limited, 
especially in oligotrophic systems where export of organic partic- 
ulate matter to depth is low. This is illustrated by the uncoupling 
between BP and PP observed in our study sites. In this study 
we demonstrated that N2 fixation and BP rates were enhanced 
when DFAA were added (Table 3). In addition to labile organic 
compounds such as DFAA, the mesopelagic aphotic layer is 
rich in organic carbon in the form of planktonic microgelati- 
nous marine snow (Azam, 1998; Azam and Malfatti, 2007). This 
includes protobiofilms that are transparent exopolymer parti- 
cles (TEP), heavily colonized by bacteria (Bar-Zeev et al., 2011, 
2012). In oligotrophic systems these protobiofilm shuttles con- 
tribute significantly to the organic carbon pool especially at 
depths where other C sources are scarce (Azam and Malfatti, 
2007; Ward et al, 2009; Bar-Zeev et al., 2012). While these par- 
ticles are rich in C, the high C:N ratios of TEP (Engel and 
Passow, 2001; Passow, 2002) make them poor N sources and 
would seem to require the use of some N supplement to foster 
microbial growth. Heterotrophic N2 fixation could provide an 
excellent source of N to fuel microbial growth on, within, and 
near these particles. Indeed, our data demonstrate a positive rela- 
tionship between TEP concentration (.R 2 = 0.49, P = 0.07, « = 
15), TEP size (R 2 = 0.90, P = 3 x 10~ 5 , n = 15), and N 2 fixation 
(Figure 3). The ability to fix N2 could provide diazotrophic 



Frontiers in Microbiology | Aquatic Microbiology 



August 2013 | Volume 4 | Article 227 | 8 



Rahav et al. 



Marine aphotic oxygenic dinitrogen fixation 



heterotrophs with a competitive advantage over non-diazotrophic 
heterotrophic bacteria in carbon-rich and nitrogen-poor micro- 
environments such as these. Moreover, in aggregates such as 
marine snow, TEP, and other marine microgels, bacterial col- 
onization causes high metabolic consumption of oxygen and 
decreasing oxygen gradients from the surface to the anaerobic 
center (Paerl and Prufert, 1987). Reduced oxygen concentrations 
would theoretically benefit diazotrophic bacteria because nitro- 
genase is an oxygen-sensitive protein (Postgate, 1998). We spec- 
ulate these anaerobic, high C, low N microenvironments could 
promote the competitive success and activity of heterotrophic 
diazotrophs at depth. 

Our novel data thus provides a conservative estimate of an 
additional 37 to 75% source of N via aphotic N2 fixation in 
the oligotrophic Levantine basin and the GA. Our measure- 
ments substantiate early indications (Paerl and Prufert, 1987) 
and recent phylogenetic characterizations (Riemann et al., 2010) 



highlighting a pivotal ecological role for diazotrophs in deep, oxy- 
genated, N03-rich waters. Because these areas comprise the bulk 
of the world's oceans they could represent a hitherto unaccounted 
for source of new N that may contribute to global productiv- 
ity and could reduce the existing discrepancies [along with other 
newly described contributions of N2 fixation (reviewed in Knapp, 
2012)] between sources and losses of N found in global marine 
N budgets. 
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